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Abstract 

The deterministic channel model for wireless relay networks proposed by Avestimehr, Diggavi and 
Tse '07 has captured the broadcast and inference nature of wireless communications and has been widely 
used in approximating the capacity of wireless relay networks. The authors generalized the max-flow 
min-cut theorem to the linear deterministic wireless relay networks and characterized the unicast capacity 
of such deterministic network as the minimum rank of all the binary adjacency matrices describing 
source-destination cuts whose number grows exponentially with the size of the network. In this paper, we 
developed a fast algorithm for finding the unicast capacity of a linear deterministic wireless relay network 
by finding the maximum number of linearly independent paths using the idea of path augmentation. We 
developed a modified depth-first search algorithm tailored for the linear deterministic relay networks 
for finding linearly independent paths whose total number proved to equal the unicast capacity of the 
underlying network. The result of our algorithm suggests a capacity-achieving transmission strategy 
with one-bit length linear encoding at the relay nodes in the concerned linear deterministic wireless 
relay network. The correctness of our algorithm for universal cases is given by our proof in the paper. 
Moreover, our algorithm has a computational complexity bounded by 0(|Vx| ■ +d- \Vx \ • C^) which 
shows a significant improvement over the previous results for solving the same problem by Amaudruz 
and Fragouli (whose complexity is bounded by 0{M ■ \£\ ■ C^) with M > d and \£\ > \Vx\) and by 
Yazdi and Savari (whose complexity is bounded by 0{L^ ■ AI^^ ■ hf) + L ■ ■ C ■ Hq) with ho > C). 

I. Introduction 

The complex signal interactions in wireless relay networks challenge the study of wireless 
information flow for many years. To characterize the capacity and capacity-achieving trans- 
mission schemes for wireless relay networks still remain open questions. Towards this end, 
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the deterministic channel model for wireless relay networks proposed by Avestimehr, Diggavi 
and Tse [1] has been a significant progress. The broadcast and inference are two fundamental 
features of wireless communications. The deterministic channel model captures the broadcast 
and inference features of wireless communications in addition to converting the wireless relay 
networks into deterministic networks. Studying the information flow in the deterministic networks 
provides a way to find out the approximated capacity and corresponding transmission strategies 
for original wireless relay networks. 

Gaussian channel has been the most widely used channel model for the link channels in 
wireless relay networks. The deterministic channel model quantizes the transmitted signal into 
different bit levels and at the receiver keeps the signal bit levels above the noise level (which 
depends on the signal to noise ratio (SNR) of the channel) so as to convert the original Gaussian 
channel into a deterministic channel without random noise variables. The broadcasting of signal 
at the transmitter is still preserved in the deterministic channel and the interference of signal at 
the receiver is modeled by modulo two sum of the bits arrived at the same signal level. 

Now we introduce the deterministic channel model by using the example of a point-to-point 
AWGN channel from [1]. Consider an AWGN channel y = hx + z where z ~ A/'(0, 1) (A/" 
represents Gaussian distribution) and h = ySNR- Assume x and z are real numbers, then we 
can write y ^ 2" Ym=i 3^(02"* + ^'^iix{i + n) + z{i))2~^ where n = [| log SNR] . If we think 
of the transmitted signal x as a sequence of bits at different signal levels, then the deterministic 
channel model truncates x and passes only its bits above noise level. Fig. \T\ gives a concrete 
example. At the transmitter node and receiver node R^, each small cycle represents a signal 
level. Assume n = 4, so only the first four most significant signal levels or bits from are 
received at Rx- Accordingly each edge in the model can transmit one-bit information at a time. 
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Fig. 1. An example of the deterministic channel model for a point-to-point Gaussian channel. 

The deterministic channel model we discussed above is called linear finite-field deterministic 
channel model in [2], which is referred to as linear deterministic channel model in this paper. 
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In [1][2], the unicast (i.e., with one source S and one destination D) capacity C of the linear 
deterministic wireless relay networks was characterized as the minimum rank of all the binary 
adjacency matrices describing S—D cuts. Since the total number of such cuts grows exponentially 
with the size of the network, an exhaustive search for the minimum rank of the adjacency matrix 
for these cuts results in an algorithm with complexity exponential in the size of the network. A 
polynomial-time algorithm is desirable for finding the unicast capacity of the linear deterministic 
wireless relay networks. We will see later that the path augmentation algorithms for graphs cannot 
be directly applied here because the definitions for the cut value for two cases are different (see 
the definition for the cut value in a linear deterministic wireless relay network in Section UIl). 
In this paper, we only consider linear deterministic wireless relay networks. Since an arbitrary 
deterministic wireless relay network can be unfolded over time to create a layered deterministic 
network through time-expansion technique [2], our algorithm developed in this paper for layered 
networks also applies to general networks. 

In [3], Amaudruz and Fragouli proposed a polynomial-time algorithm for finding the unicast 
capacity of a linear deterministic wireless relay network by trying to identify the maximum 
number of linearly independent paths in the network using the idea of path augmentation. In 
[4], Yazdi and Savari developed a two-dimensional Rado-Hall transversal theorem for block 
matrices and used the submodularity of the capacity of a cut to formulate the problem as a 
linear program over the intersection of two polymatroids so as to solve the problem in polynomial 
time. Compared with these previous results, our algorithm has a significantly less computational 
complexity as explained in Section IIV-AI Moreover, our algorithm comes with a more intuitive 
understanding as we explained in Section Hill 

The paper is organized as follows. In Section [III we introduce the notations and definitions 
used throughout this paper. Section HlI] gives a detailed description of our algorithm for finding 
the unicast capacity of any given linear deterministic wireless relay network. Section |IV] is about 
the algorithm analysis including the proof of correctness and complexity analysis. Section |V] 
concludes the paper. 

II. Notations and Definitions 

Let Q = {J\f, V, S) be a layered deterministic wireless relay network [2] where J\f is the set of 
super nodes referring to the nodes in the original wireless relay network, V is the set of nodes 
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referring to different signal levels incident to super nodes and E is the set of directed edges going 
from node to node. Denote V = U Vy with = {x : (x, y) E £} and Vy = {y : (x, y) E £}. 
We shall call Vx as transmitting nodes and Vy as receiving nodes. For example, in Fig. [H 
J\f = {Tj., Rx} and has five transmitting nodes and has five receiving nodes. In a layered 
network, all paths from the source to the destination have equal lengths [2], so we can divide 
the set of super nodes into different layers each layer containing only super nodes with the same 
distance to the source. Assume there are L layers of super nodes in Q and M is the maximum 
number of super nodes in each layer. The source super node S consists the first layer and the 
destination super node D consists the last layer. Let ^/{xi) (or J\f{yj)) denote the super node 
where a transmitting node Xj (or a receiving node yj) belongs to. Let C{N) (or C{xi), C{yj)) 
denote the layer number where super node (or Xi, yj) belongs to. In a layered network, if 
{x,y) E £, {x',y') E £ and C{x) = C{x'), then we must have C{y) = C{y') = C{x) + 1 = 
Cix') + 1. 

A cut, n, in a deterministic relay network ^ is a partition of the super nodes J\f (together 
with their incident nodes) into two disjoint sets f2 and 0,^ such that S E Cl and D E ^l''. For 
convenience, we call a cut a layer cut if all edges across the cut are emanating from nodes 
belonging to the same layer, otherwise we call it a cross-layer cut. We say an edge (xj, yj) E £ 
belongs to layer cut / if C{xi) = I. For a layered deterministic network, there are exactly L — 1 
layer cuts. 

The adjacency matrix T for a set of transmitting nodes x = {xi,X2, .■■Xm},Xi E Vx, and a 
set of receiving nodes y = {yi,y2, ...yn},yi G Vy in a deterministic relay network is a binary 
matrix of size m x n with rows corresponding to {xi,Xi E x} and columns corresponding to 
{Vi^yi £ y} and T{i,j) = 1 if {xi,yj) E £. The adjacency matrix for a set of k edges is the 
binary adjacency matrix for the set of their transmitting nodes and the set of their receiving 
nodes. 

A set of k edges are said to be linearly independent if the adjacency matrix for them has rank 
k, otherwise they are said to be linearly dependent. In a layered deterministic network, each 
S — D path is of length L — 1 and crosses each layer cut exactly once. A set of A; 5" — D paths 
are said to be linearly independent if each set of their edges of size k crossing each layer cut 
are linearly independent, otherwise they are said to be linearly dependent. Lemma |4] shows that 
k linearly independent S — D paths correspond to a transmission scheme of rate k. 
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In a deterministic wireless relay network, there are intermediate super nodes (exclude S and 
D) corresponding to the relay nodes in the original wireless relay network which have both 
transmitting nodes and receiving nodes. It is shown in Lemma |4] and Theorem [2l that one-bit 
length linear encoding functions at the relay super nodes are sufficient for constructing capacity- 
achieving transmission schemes between S and D for the underlying linear deterministic relay 
network. 

Let Sfi be the set of edges crossing the cut Q in a linear deterministic relay network. The 
cut value of Q in the linear deterministic relay network is defined as the rank of the binary 
adjacency matrix for Sn, which equals the number of linearly independent edges in Sq. Note 
that here the cut value defined for linear deterministic wireless relay networks is different from 
that for graphs (which is just the number of edges crossing the cut). It is proved [1][2] that the 
unicast capacity of a linear deterministic wireless relay network is equal to the minimum cut 
value among all cuts. 

III. Our Algorithm 

A. Algorithm Outline 

The max-flow min-cut theorem has been generalized to the linear deterministic wireless relay 
networks, but since the definitions of cut value for graphs (the number of edges crossing the cut) 
and for linear deterministic relay networks (the number of linearly independent edges crossing 
the cut) are different, the path augmentation algorithms for the latter are different from that for 
the former. As we will see later from this paper, the similarity is that in both cases the path 
augmentation algorithms operate in iterations and in each iteration they complete an additional 
path. The difference is that contrast to the simple addition of an available edge to a path in 
the path augmentation algorithms for graphs, the addition of edges to a path in a deterministic 
network has to satisfy some rank requirement of the adjacency matrix to avoid linear dependence 
among used path edges in each layer cut. 

Our algorithm is basically a path augmentation algorithm for finding the maximum number 
K of linearly independent S — D paths in a layered linear deterministic relay network, where 
K proves to be equal to C in Section IIV-B[ Moreover, the K identified paths correspond to a 
capacity-achieving transmission strategy for the underlying deterministic network. 
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The algorithm operates in iterations. In iteration k, a modified depth-first search (MDFS) 
algorithm tailored for the linear deterministic relay networks is carried out on the graph Q trying 
to find an S — D path in addition to the k — 1 paths found in the first k — 1 iterations (stored in a 
structure V) so that the k found paths (stored in a structure V) are linearly independent. If MDFS 
returns True indicating that the A;th S — D path is found, then totally k linearly independent paths 
have been found and the algorithm continues to iteration A; + 1 . If it returns False indicating that 
no S — D path is found, then no more independent path exists and the algorithm stops while 
those identified paths in V suggest a capacity-achieving transmission scheme. 

B. Preliminaries 

The following notations or structures are used in our algorithm. V and V' are structures storing 
information about the k — \ paths found in previous k — \ iterations and information about the 
updated k — 1 paths and the partial path found in the current fcth iteration respectively. Denote 8u 
(£li) as the set of edges in £ used by paths in V (in layer cut i). The MDFS algorithm ensures that 
edges in S^^ are linearly independent, i.e., rank(T(£^^)) = \£l\. Denote V^.^ = {x : {x,y) G £1} 
and = {y : {x,y) G S^}. £u (SI), V*„, and V are initialized at the beginning of every 
iteration according to V and maintained by MDFS in the current iteration subject to changes. 
As we will notice in Section |III-C[ each move MDFS makes ensures that rank(T(£^*)) = \£l\. 
Let £1 ^ , V* , ^ and V!,„ ^ denote the instantaneous sets of £1, V!„ and Vi, when a transmitting 
node Xj is being explored in layer cut i = C{xj) by MDFS in a certain iteration. Furthermore, 
denote Vx,n = {x : x G and ^{x) = N} and Vy^N = {y : y E Vy and A/'(?/) = A^}. Denote 
£-p = {e : e & £,e used by some paths in V} and Vxv = {x '■ [x, y) G £-p for some y} and 
Vy-p = {y : (a;, y) G £v for some x}. 

Definition 1: Vxlpan- We define Vxipan as the set satisfies Vxlpan C V^„^^ and 

n^.,v;„,.;= n^^lu..,) (1) 

X A 

xspan 

Since T(V*„ ^.^ , V*„^^ ) has full rank, Vxipan is the unique such set. 

Let function Span(xj) be the function for computing Vxipan- In iteration k + 1 of our algorithm, 

IV* I = k. 

I XU,Xj I 

Lemma 1: Let |V^„^. | = k = \Vy^^.\. The computational complexity of Span(xj) for finding 
the set Vxipan is bounded by Oi^). For Vx G Vxipan, rank(r(V^„_^^,, V^„^^^)) = rank(T(V;„^^^, + 
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^j-x,Hu,.,)) = A;andrank(T(VlVn+a:,,V;,,,^)) = rank(T(VSpan, V;„,,^,)) = rank(T(V:ipa,n+ 

Xj — X,Vy^;^^)) = \Vxspan\ < k. 

Proof: To solve the set Vxipan is equivalent to solving the system of linear equations, 
Vxipan ■ T{V:u,x,,Ku,x,) = T{xj^Ku,x,) GF{2) which cau be accomplished in time Oi^) 
by using Gaussian elimination. The second statement is obvious. ■ 
Consider the subgraph consisting of nodes xj U V*.„ ^.^ U V*„ ^.^ and the edges connecting them 
in Q. Let ^^^^ denote the graph obtained by reversing the directions of the edges in E^^^^. in the 
above subgraph. 

Lemma 2: There is a directed path from xj to any x G Vxspan in graph Qg^i,. Let FindlndPaths 
be the function for finding out all these |Vxipan| paths from xj. The computational complexity 
of FindIndPaths(xj, Vxipan) is bounded by 0{k^) in iteration k. 

Proof: From Lemma [H for Vx G Vxipan, rank(r(VL,x, , ,^^ )) = rank(r(VL,x, + - 
X, V*„^^)) = k where k = \V\. Introduce an auxiliary receiving node y' and an edge (x, y'). It's 
easy to see that rank(T(V;„,,^, + x,-, V^„,,^, + y')) = k + 1. Given rank(r(V:,,,^,, V^„,,^,)) = k 
and rank(T(V*.„^^ + Xj,Vy^ .^. + y')) = k + 1, there is a size k prefect matching between V*„^^, 
and V*„2,^, Ml = f"^^^ being such a matching, and a size A; + 1 perfect matching between 
'^xu,xj + ^^'^ '^yu,xj + y'- Using a similar argument as in finding the maximum bipartite 
matching, we claim that there is an alternating path, relative to the matching Mi, starting from 
an unused transmitting node Xj to an unused receiving node y', alternating between edges 
not in the current matching Mi and edges in the current matching Mi, i.e., there is a path 
Px.^y' = {{xj,yi),{yi,xi),{xi,y2),{y2,X2),...{xm-i,ym),iym,Xm),{xm,y') = ix,y')} with 
< i < m being edges in Mi = S^^^^. So we proved that there is a path Px^^x = 
{ixj,yi),{yi,xi),{xi,y2),iy2,X2),...{x^-uym),iym,Xm) = iym,x)} with {xi,yi),l < i < m 
being edges in i^*^^. Without loss of generality, we also use the following representation of 
the path P^^^x = {{xj,yi), {xi,yi), {xi,y2), {x2,y2), ■■■{xm-i,ym), {xm,ym) = {x,ym)} with 
(xj, 1 < ? < m being edges in S^,,... 

In iteration k of our algorithm, |V*,„| = |V*„| = A; — 1, so the number of nodes in Q^,^^ is 
bounded by 0{k), which also means the number of edges in Q^^^ is bounded by 0{k'^). To find 
directed paths from Xj to all x G Vlipan in takes time bounded by 0{k'^) by using some 
well-known graph traversal algorithms, like breadth-first search. ■ 
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Let's briefly recaU the depth-first search (DFS) algorithm first. DFS algorithm is a well-known 
algorithm for traversing graphs. To find a path, DFS progresses by exploring the outgoing edges 
of a node before exploring any other outgoing edges of the node's predecessor. In this way, it 
proceeds deeper and deeper until it reaches the goal node when it stops or until it reaches a node 
without any outgoing edges when it backtracks to the most recent node that it hasn't finished 
exploring. In the search process, each node in the graph would be explored at most once. 

Our modified depth-first search (MDFS) algorithm for linear deterministic relay networks 
inherits the basic forwarding and backtracking operations from the basic DFS algorithm and 
each super node in A/" is treated like a node in DFS. The difference is that in order to find 
an S — D path in the graph Q = {J\f, V,£) linearly independent to the paths in V, we have 
to avoid linear dependency between different paths in our algorithm, which means that instead 
of allowing a valid forwarding move along each outgoing edge when a node is explored as in 
DFS, more constraints should be imposed on the forwarding moves of a super node in MDFS. 
We propose the following MDFS algorithm to accomplish the path augmentation task in our 
problem. 

C. Modified Depth-First Search (MDFS) Algorithm 

In the following, the exploration to a super node N E J\f or a node i; G V by MDFS refers 
to that MDFS has extended the path found in the current iteration to N or J\f(v) and now it 
continues to extend the path further from N or from v. The exploration to a super node is 
realized as the exploration to some of its unexplored incident nodes v (which we call admissible 
nodes of in Definition |2l). Once a super node or a node is explored, it's labeled explored. 
During the running time of MDFS, each super node N E M could be labeled as unexplored 
or explored indicating that it hasn't or has been explored by MDFS. Each node in V could 
be labeled as unexplored or explored indicating that it doesn't allow or allows exploration by 
MDFS. There is also a type labeling with each node in Vx indicating how it can be explored 
by MDFS if it allows exploration. Let SetLabel(X, LABEL) be the function setting label of 
X (X can be a super node or a node) as LABEL (LABEL can be explored or unexplored). 
Let LABEL=GetLabel(X) be the function returning the label of X. Let SetType(x, B) be the 
function setting type of transmitting node x as B (B can be 1,2, or 3). Let B =GetType(x) be 
the function returning the type of x. 
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Definition 2: Admissible nodes and admissible forwarding moves of super nodes - We 

define the admissible nodes for a super node when is explored by MDFS as follows: 
V^'^ = {x : x e Vx,N,x ^ Vxii'^\x labeled unexplored or y : y e Vy^NiV ^ ^yViV labeled 
unexplored}. The exploration to a super node N is realized as the exploration to its admissible 
nodes. The following forwarding moves starting from V^'^ are defined as admissible forwarding 
moves allowed in MDFS algorithm on the graph Q = (A/", V, £) when it explores a super node 
N. 

1) Type 1: moving forward along {x,y) E £ with x E V^, GetType(x) = 1,2, or 3 and 
y ^ V^i'^^ rank(T(vfi? + x, v£,x^ + y)) = \V\ + 1 and GetLabel(Ar(2/)) =unexplored. 

2) Type 2: moving forward along the path Px^x' (as proved to exist in Lemma [2l) for any 

e ^xspan with X E V^^ , GetType(x) = 1 or 3. 

3) Type 3: moving forward along {x,y) E 8 with y E V^, {x,y) E Su}x^^~^. 

When is explored, x E Vj^ with GetType(x) = 1 or 3 would allow to start type 1 and type 
2 admissible forwarding moves, x E with GetType(a:) = 2 would allow to start type 1 
admissible forwarding moves only and y E would allow to start type 3 admissible forward 
moves. There is an ordering in exploring V^'^: type 2 and type 3 nodes should be explored first, 
then type 1 nodes and finally the receiving nodes. 

Definition 3: Modified depth-first search algorithm (MDFS): The MDFS algorithm is 
defined in terms of initialization, exploring of a super node, labeling of A/", V and updating 
of V and £u as follows: 

1) Initialization: Set J\f and V as unexplored and Vx as type 1 nodes. 

2) Exploring of a super node N: When MDFS explores a super node A^, it means that a 
partial path P\v\+i from to A^ has been found in addition to the |P| complete S — D 
paths. As mentioned in Definition [21 the exploration to a super node A^ is realized as 
exploration to its admissible nodes V^'^ and three kinds of admissible forwarding moves 
from V^'' are allowed. There is an ordering in exploring V^'': type 2 and type 3 nodes 
should be explored first, then type 1 nodes and finally the receiving nodes. Once a super 
node or a node is explored by MDFS, it is labeled explored. Now we explain how to 
understand these three admissible forwarding moves and how to label J\f, V and update 
V, £u with these moves. 
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a) Type 1 admissible forwarding move: A type 1 admissible forwarding move can be 
understood as that MDFS extends P\v\+i along an edge (x, y) from super node 

to super node Af{y), i.e., P\v\+i = P\v\+i + i^.y). V is updated accordingly. Then 
MDFS makes M^y) with C{J\f{y)) = C{N) + 1 the next super node to be explored 
when Su^'^^ is updated as Su^^'' + {x,y) and \Su^^''\ increases from \V\ to \V\ + 1. 

b) Type 2 admissible forwarding move: A type 2 admissible forwarding move can be 
understood as that MDFS updates V according to the path Px^x' (proved to exist in 
Lemma [2l) for some x' G Vxspan follows. Let Px^x' = {{x,yi), {xi,yi), {xi,y2), 
■ ■■{xm, ym) = {x', ym)} with {xi, yi) , 1 < i < m. being path edges used hy m S - D 
paths in V existing before the current move, denoted as Pi,l < i < m. After the 
current type 2 forwarding move, Pi,l < i < m and P\v\+i are updated to P- and 
P'-p^^i as follows. Let Pi(Ni, N2) denote the segment of Pi from super node A^i 
to N2. Pi = Piv\+i{S,^{x)) + (x,yi) + Pi{^{yi),D), PI = Ar(x,_i)) + 
(xi^uyi) + PiiAfiy,),D),l <% <m and P[^^^^ = P„(5, Ar(x')). V is updated 
accordingly. After the type 2 admissible forwarding move, x' is labeled as unexplored 
and set as type 2 node. Then MDFS makes Af{x') with C{Af{x')) = C{N) the next 
super node to be explored when Su^^'' is updated as Su^^'' + {x,yi) — {xi,yi) + 
{xi,y2) - {x2,y2) + ••■ + {xm.~i,ym) " {x' , ym) and l^^^-*! keeps to be \V\. 

c) Type 3 admissible forwarding move: A type 3 admissible forwarding move can be 
understood as that MDFS updates V' along the edge (x, y) as follows. Let Px be the 
path in V existing before the current move where {x, y) belongs. If Px 7^ P\v\+i^ 
then after the current move, MDFS updates Px and P\v\+i to and P\-p\j^i as 
follows. P'^ = P\v\+i{S, N) + Px{N, D) and P,'^,^, = P,(S, A/'(x)). If Px = P\v\+i, 
then after the current move, MDFS updates P\v\+i to P'-p^^^i as follows. P|'p\_^_l = 
P\-p\+i{S , Af (x)) . V is updated accordingly. After the type 3 move, y is labeled 
explored, x is labeled as unexplored and set as type 3 node. Then MDFS makes 
J\f{x) with C{N'{x)) = C{N) — 1 the next super node to be explored when 8u^'^^~^ 
is updated as £u^^^~^ — {x,y) and \Su^^''~^\ decreases from \V\ + 1 to \V\. 

Lemma 3: In iteration A; + 1 of our algorithm, MDFS defined in Definition [3] maintains a set 
of size k linearly independent S — D paths while it tries to complete a(A; + l)thS' — D path 
Pk+i- And when MDFS explores a super node A^, it means MDFS extends P^+i from 5 to 
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in addition to k complete S — D paths (all stored in V') and there are A; + 1 linearly independent 
used path edges in each of the first C{N) — 1 layer cuts and k linearly independent used path 
edges in each of the rest layer cuts. These two facts lead to the conclusion that when MDFS 
returns True in iteration k + 1, totally k + 1 linearly independent S — D paths are found (and 
stored in V'). 

Proof: We prove by induction. Assume that in the first k iterations of our algorithm, MDFS 
already finds k linearly independent S — D paths. Then at the beginning of iteration k+1 before 
we call MDFS, these two statements are clearly true. Now it's sufficient for us to show that 
these three kinds of forwarding moves allowed in MDFS keep these two statements true after 
each forwarding move of MDFS. First, clearly a type 1 forwarding move along an edge (x, y) 
increases |£^n''^''| from k io k + 1 with rank(T(£^u *''^^)) = A; + 1 and it has no impact to the set 
of k S — D paths existing before the current move. So these two statements remain true after 
a type 1 forwarding move. Second, a type 2 forwarding move keeps unchanged with 

rank(T(£^^*^^'')) = k (according to Lemma [T]) while it updates the paths Pi,l < i < m and Pk+i 
to P/, 1 < z < m and Clearly the rest k — m S — D paths existing before the current move 
are not impacted by the move. So these two statements remain true after a type 2 forwarding 
move. Third, a type 3 forwarding move extends P^+i along an edge (x, y) from super node 
Af{y) to Af{x). Because MDFS reaches Af{y) before the current move, {Su'^^'^'^'^l = k + 1 and 
after the current move, ^''^''^■'l = k with {x,y) deleted. Clearly MDFS maintains k linearly 
independent S — D paths after the current move when it proceeds to explore J\f{x). So these 
two statements remain true after a type 3 forwarding move. ■ 
Lemma 4: Let V be the set of S — D paths found by our algorithm and \V\ = K, then the 
paths in V correspond to some transmission scheme of rate K from S to D for the underlying 
deterministic network. 

Proof: We prove by constructing a transmission scheme of rate K from S to D for the 
underlying deterministic network by using the K paths in V. Lemma |3] says that these K paths 
in V found by our algorithm are linearly independent. Let V^^y C Vy^N (V^at ^ yx,N) be the 
set of transmitting (receiving) nodes incident to N that are used by some paths in V. Clearly 
|V^7vl = l^x^Afl- The relay function at each relay super node N could be any one-one mapping 
between V^j^ and V^jy i-^-' ^^ch received bit from V^^y is transmitted forward by a unique 
transmitting node in Vj"^ specified by the mapping. Clearly each such mapping corresponds to 
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a full rank adjacency matrix. For simplicity, we can treat those relay functions as intra-layer 
paths which together with the paths in V (treated as inter-layer paths) completely specify a 
transmission scheme from S io D. Since in each of the L — 1 inter-layers and each of the L — 2 
intra-layers the adjacency matrix for used path edges has full rank K, the overall transfer matrix 
from S to D along these paths is the product of all these adjacency matrix and also has full 
rank K which guarantees that transmission information rate K is allowed between S — D with 
the transmission scheme defined by these K paths in V and any one-one mapping function for 
each relay super node. ■ 
Let's justify the rule that type 2 nodes shouldn't start any type 2 forwarding moves in MDFS. 
By definition, a type 1 or type 3 transmitting node x E can start type 1 or type 2 admissible 
forwarding moves of while a type 2 transmitting node x' E can only start type 1 admissible 
forward moves of A^. Lemma [H helps in explaining why type 2 transmitting node shouldn't start 
any type 2 forwarding moves. When x' is labeled as type 2 node after a type 2 forwarding 
move along the path Px^x' for some type 1 or type 3 node x, it means that x' E Vf^^^^. It's 
easy to see that when x' is explored right after the type 2 forwarding move along path Px_»zs 
'^x'span = '^xspan + X — x' . Sincc X is cxplorcd and x is of type 1 or type 3, all transmitting 
nodes in V^^p^^ could be explored following type 2 forwarding moves from x. It means that all 
transmitting nodes in V^^p^„ are guaranteed to have the chance of being explored by MDFS as 
long as X is explored and all type 2 forwarding moves from x are allowed. So type 2 node x' 
shouldn't start any type 2 forwarding moves to avoid redundancy. 

D. Implementation Details of MDFS 

The MDFS algorithm is implemented in the (Res,P') =MDVS{g,V,V' ,N,D,k) function. 
The input parameters are the graph Q, information about the k—\ found S — D paths stored in 
structure V, information about the updated k — 1 S — D paths and partial path Pk stored in a 
structure V , starting super node from where to complete Pk, ending super node D where to 
end Pk and the iteration number k. In each iteration, we call MDFS by initializing to be 5" for 
finding an S* — -D path. The function returns (True, V) with \V'\ = \V\ + 1 = k if the A;th S — D 
path is found in current iteration and (False, V) with V' = V li no path is found. Our algorithm 
calls MDFS(^, V, V, S, D, k) in iteration k for increasing k until MDFS(^;, V, V, S, D, K + 1) 
returns Res=False for some K, then our algorithm stops and claims that K = C. Each time 
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MD¥S(Q,V,V' , S, D, k) returns Res=True, V is updated as V and used in next iteration. 



E. Accelerating The MDFS Algorithm 

From definition, we see that the main computational complexity of the MDFS algorithm comes 
from the rank computation of the associated binary adjacency matrix in deciding whether an 
edge (x, y) qualifies for a type 1 admissible forwarding move when x of type 1, type 2 or 
type 3 is explored and from the computation of V^gpan deciding whether a type 2 admissible 
forwarding move exists between a pair x, x' when x of type 1 or type 3 is explored by MDFS. In 
this subsection we explore some useful combinatorial properties related to the rank of the binary 
adjacency matrix and develop an equivalent but computationally simple method to replace the 
rank computation in MDFS. 

In iteration + 1 of our algorithm, let Y^^"^ = {y ■ C{y) = C{xj) + 1 = i + l and y ^ ^.^ }, 
let yj+i = {y:ye Y^+^ and (x, y) e S for some x e {xj U VSpan}} and let Y^+j\ = {y y e 
F^+i and rank(T(VL,,^, + x„ V^„,,^. + y)) = k + 1}. 

Lemma 5: Y^^j^. is a subset of Y^^}. The rank check of the adjacency matrix (rank(T(V*„^^ + 
'^yu.xj+y)) = kov k+1) for V?/ G is equivalent to checking T{xj, y) = Ex-gv,%,„ ^(^' 
or not, whose computational complexity is bounded by 0{k). 

Proof: First we prove that for G Yj:^:^^, if there is no (x, y) E S with x G {xj U Vxipan}, 
then T(V*.„ + Xj, Vl^ ^^. +y) = k. If there is no (x, y) E S with x G {xj U Vxipan}, then the 
vector r(V*.„^^ + x^ , y) has zero entries on rows indexed by Xj U Vxipan, together with equation 
O we have T{x„ V;„,,^ +?/) = E..ev:.V„ ^(^' Hu,x,+y)- So rank(r(V;,,,^. + x„ V^^^^^+y)) = 
rank(r(V;„,,^, V;,,,^ + y)) = k. So we proved ^4+/, is a subset of y^+i. 

Weknowrank(T(V:,,,^., V^„,,^+y)) = and |V;,,,J = /c. If rank(T(V:,,,^.+x„ V^„,,^.+2/)) = 

x' 

k, then it must be T(xj, V^„^^^. + y) = E ^ T{x,Vl^^^^ + ^) ^^'^p"" ^ 

Given equation ([B and VSpan is unique, we have Vxspan = Vxipan and T(xj , V*„^^ + ?/) = 
S ^(^' + equivalent to T(x„ y) = E.^v^.^ ^(^' So rank(r(V;„,,^ + 

Xj, Vi^,x, +y)) = k or k + l is equivalent to T{xj, y) = E^gv^.V^ ^(^' °^ " 
Lemma [5] says that given Vf^p^^ we can simplify the rank computation of the associated binary 
adjacency matrix in deciding whether an edge (x, y) qualifies for a type 1 admissible forwarding 
move when x of type 1, type 2 or type 3 is explored. We also know that when x' is explored 
as a type 2 node right after the type 2 forwarding move along path Px^x' for some x of type 1 
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or type 3, V^^^^ = V^,p„„ + x-x', so we don't need to compute V^'^p„„ for any x' of type 2 
being explored as long as Vf^p^^ is computed. 

Lemma 6: In the last iteration K + 1 of our algorithm when no more S — D path is found 
by MDFS, all Af{xf,) with x^ E Vxipan and all ^{y) with y G Y^.^fl must have been explored 
by MDFS before MDFS returns if Xj has ever been explored. 

Proof: Since there is no more S — D path exists in the last iteration of our algorithm, all 
possible admissible moves would be tried by MDFS. It means each A/'(xfc) with Xk G VSpan 
would be explored by following a type 2 forwarding move along the path Pxj^Xk if Xj of type 1 
or type 3 has ever been explored. For any Xi of type 2 being explored right after a type 2 move 
along the path Px^^xi starting from some Xj of type 1 or type 3, we know that Xj is explored 
and V^ipan = Vxipan + Xj - Xi. So if all Af{xk) with Xk G Vxspan are explored, all J^{x'f.) with 
x'k e VSpan are also explored. 

From Lemma [51 for \/y G Y^.jl, there is (x, y) G £^ for some x G {xj U Vxipan}- We already 
proved that all ^{xk) with Xk G Vxipan would be explored by MDFS in the last iteration of our 
algorithm if Xj has ever been explored. So the edge (x, y) with x G {xjVJVxipan} and y G X^f^j^^ 
must be considered by MDFS in the last iteration given that Xj has ever been explored. If x = Xj, 
(x, y) allows a type 1 forwarding move when Xj is explored and then ^{y) will be explored. 
Now assume x G Vxipan- Given Xj is explored, x is also explored and V*„ 3, + x = V*„^^, + Xj 
and Vi^^^ = Vl^^^. Given y G yjj,^ i.e., rank(r(VL^^^. + Xj^V^^^^^ + y)) = + 1, we have 
rank(T(V*„ 2. + x, V*„^ + ?/)) = /\ + 1. Then (x, y) allows a type 1 forwarding move when x 
is explored and then M{y) will be explored. So all M{y) with y G niust be explored by 

MDFS in the last iteration of our algorithm if Xj has ever been explored. ■ 

Table HI gives a pseudo-code description of the MDFS algorithm described above. Refer to [5] 
for an implementation of the algorithm proposed in the current paper. 

IV. Algorithm Analysis 

A. Complexity Analysis 

Theorem 1: MDFS algorithm defined in Section UlI] terminates in finite time as the total 
number of explorations to transmitting nodes invoked by MDFS is bounded by OdV^:! ■ k) 
in iteration k of our algorithm. The total computational complexity of our algorithm is bounded 
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TABLE I 



(Res={Trae, False},P')=MDFS(e, V, V , N, D, k) 
SetLabel(A'^, explored) 

for G Vx,N with x Exu^'' and GetLabel(a-) = unexplored and GetType(x) = 2 
SetLabel(a::, explored) 

for Ve = {x,y) £ C/.incidentEdgeTypel (x) and GetLabel(A/'(i/)) = unexplored 



' if r(a.,y)/E.'6v,x^^_n^',J/) 



ff(~'-ff''^'+e; Update(P') 



if M{y) = D { return (True,P') } 
{ J (Res,P') = MDFS(e,P,P',A/'(y),-D,fc) 



else < 



I if Res = True { return (True.'P') } 



[ [ [ - ff*'^' - e; Restore(P') 



for Va; £ Vx,n with x Exu^^ and GetLabel(a;) = unexplored and GetType(x) 7^ 2 

(first for X with GetType(a;) = 3, then for x with GetType(x-) = 1) 
f 

SetLabel(a::, explored) 
^ispan = Span(x) 

for Ve = (a-, G C/.incidentEdgeTypel (a;) and GetLabel(A/'(y)) = unexplored 



[ [ - - e; Restore(P') 

■P^^v-^p^,. = FindIndPaths(a;, V^3p^„) 

for Mx' G V^spa„ with Px-,x' = {ei, 62, ...62™} = {(s, yi), (yi, xi), (si, ^2), •••(y,n, a;™) = (j/„i,a:')} 
SetLabel(a:', unexplored); SetType(x', 2); v£p„„ = V^sp„„ - x + x 

^ fif + ei - e2 + 63 - e4 + ... + e-im-x - e2m; Update(P') 
< (Res,P') =MDFS(e,P,P',A/'(a:'),i3,fc) 

if Res = True { return (True,P') } 
[ Eu^^^ ^ Eu^^'' - ei + 62 - 63 + 64 - ... - 62m-i + 62m; Restore(P') 
for Ve = (a;, y) G C/.incidentEdgeType3(Af) with GetLabel(2/) = unexplored 
SetLabel(y, explored) 
SetLabel (a;, unexplored); SetType(a;, 3) 




cC.{N)-\ ^£{iV)-l , ^ 

E-a ^ Eu — e; Update^P ) 



(Res,P') = MDFS{g,r,r',Af{x),D,k) 
if Res — True { return (True.'P') } 



^ £f'^'-'-£f<^'-Ve; Restore(P') 



return (False,P') 



15 



by 0(|Va;| ■ + d - \ Vx\ ■ C^) if our algorithm stops after finding C linearly independent S — D 
paths. 

Proof: In MDFS, a transmitting node is explored or re-explored exclusively as an unex- 
plored type 1, type 2 or type 3 node. Let ki, k2 and be the total number of labeling of a 
transmitting node as unexplored type 1, type 2 and type 3 node respectively. To prove MDFS 
algorithm terminates, we can equivalently prove that ki, k2 and ks are finite. 
Claim 1: ki < |V^|. 

Proof: All nodes in Vx are labeled as unexplored type 1 nodes at the beginning of the 
iteration. After an unexplored type 1 node is explored, it's labeled as explored and never relabeled 
as unexplored type 1 node again. So A;i < |Vi:|. ■ 
Claim 2: k^ < \Vx\- 

Proof: A node is labeled as unexplored type 3 node only after a type 3 admissible forwarding 
move. We next show that the number of type 3 admissible forwarding moves is bounded by \Vx\- 
A type 3 admissible forwarding move only happens on an edge (x, y) with y E Vy-p unexplored 
and y is labeled as explored after a type 3 move and never relabeled as unexplored again. Since 
\^yv\ = \^xv\ < so the total number of type 3 moves is bounded by \Vx\- The total number 
of labeling of transmitting nodes as unexplored type 3 nodes is therefore bounded by |Va;|. ■ 

Claim 3: k2 < k ■ [ki + k^). 
Proof: A node is labeled as unexplored type 2 node only after a type 2 admissible forwarding 
move. We next show that the total number of type 2 moves is bounded hy k2 < k ■ [ki + /^s). A 
type 2 admissible forwarding move only starts from some type 1 or type 3 node being explored. 
The total number of type 1 or type 3 nodes explored by MDFS is bounded by 2\Vx\- When a 
type 1 or a type 3 node x is explored, the total number of type 2 admissible forwarding moves 
that it starts is bounded by |Vf^pa„| < k. So the total number of type 2 moves is bounded by 
2 ■ k ■ \Vx\- The total number of labeling of transmitting nodes as unexplored type 2 nodes is 
therefore bounded by 2 ■ A; ■ |Vx |. ■ 

Now it is easy to conclude that the total number of explorations to transmitting nodes invoked 
by MDFS is bounded by 0(|Vx| ■ k) in iteration k of our algorithm, so MDFS terminates in 
finite time. 

Now let's consider the computational complexity of our algorithm. We already proved that 
in iteration k the total number of type 1, type 2 and type 3 transmitting nodes explored by 
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MDFS is bounded by IV^:!, 2|Va;| ■ A; and \Vx\ respectively. The worst case in computation for all 
these transmitting nodes is that (1) MDFS computes V^^p^^ (with complexity bounded by 0{k^) 
based on Lemma [T]) and finds the paths Px^v^^p^^ (with complexity bounded by 0{k'^) based 
on Lemma [2l) for any x being explored as type 1 or type 3 node and (2) MDFS checks the 
rank of the binary adjacency matrix associated with each incident edge to x (with complexity 
bounded by 0{k) based on Lemma |5]) for any x being explored as type 1, type 2 or type 3 
node. Clearly, the total number of such edges is bounded by 2(iA;|V2:| in iteration k. So the 
computational complexity of MDFS in iteration k is bounded by 0(|Vx| ■ k^ + d- \ Vx\ ■ k'^)- The 
total computational complexity of our algorithm is bounded by OdV^I -C^ + d- \ Vx\ -C^) if our 
algorithm stops after finding C linearly independent S — D paths. ■ 
Our algorithm shows a significant improvement in terms of computational complexity over the 
algorithms for solving the same problem in [3] by Amaudruz and Fragouli (whose complexity 
is bounded by 0{M ■ \S\ ■ C^)) and over the algorithm in [4] by Yazdi and Savari (whose 
complexity is bounded by 0{L^ ■ M^^ ■ + L- ■ C ■ hi)). Note that here M > d (because 
each transmitting node can transmit at most one bit information to each super node by definition 
of the deterministic channel model), \S\ > |V^| (because of broadcasting) and > C (because 
C cannot be larger than the maximum number of transmitting nodes among all layers). 

B. Proof of Correctness 

Theorem 2: Our algorithm finds C linearly independent paths in a linear layered deterministic 
relay network Q where C is the unicast capacity (or the minimum cut value among all cuts 
separating the source from the destination) of Q. 

Proof: We prove Theorem |2] by proving that when our algorithm stops the number of paths 
we find equals some cut value in Q. 

If MDFS returns Res=False in iteration A' + 1 of our algorithm, then it finds the maximum 
number K of linearly independent S — D paths in Q and we claim K = C. Consider the cut 
^Ik separating the super nodes labeled explored from the super nodes labeled unexplored in 
iteration K + 1 when the algorithm stops so that 5" G ^Ik- Clearly ^Ik is a cut separating 5" 
from D and S G Hk, D G ^2^. We prove Theorem [2] by proving that this cut value equals K, 
i.e, mnk{T{Sn^)) = K. 

Let = {{x,y) : {x,y) G E-p and C{x) = i}. Let = {x : {x,y) G and = {y : 
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{x,y) e S^}. We divide the set into four subgroups: S^^ = {{x,y) : {x,y) E S^,Af{x) E 

nK,^fiy) E fi^}, ^2 = {ix,y) ■ i^.y) e £ipM{x) e ^K^iy) e nj,}, ^3 = {ix,y) 

ix,y) E £jpM{x) E niM{y) e n'j,} and 4^ = {{x,y) : {x,y) E £'^,M{x) E fi^,Ar(2/) E 
VLk}- We divide the sets V*^ (V*p) into four subgroups accordingly, V'^j 0^lvi>^ 1 < J < 4. 
Clearly, the subgroups Vl-pj, 1 < j < 4 are disjoint, so are true for subgroups Vy-pj, 1 < j < 4. 
Denote \£'!pi\ = \ Vl'pi\ = \Vy-pi\ = Ki. Clearly Ki is the number of our identified paths (or path 
edges) that cross the cut VLk in layer cut i. Denote \S'!p2\ = Ki2, ISp^l = and \£p4^\ = K^. 
Clearly, Ki + Ki2 + Ki^ + = K. 

Let = {{x,y) : {x,y) E £,M{x) E VLK,N'{y) E VL% and C{x) = z},l < i < L, i.e., 
81^^ is the intersection of the cut VLk and layer cut i. Consider the adjacency matrix T(Efij^) 
for ^Ik- It is a block diagonal matrix with each block T^^, being the adjacency matrix for S^^^ 
and rank(T(^n,)) = Ef=i' rank(T^^). Let V^^^ = {x : {x,y) E S^J and V^^^ = {y : 
(x, y) E S'nJ. We also divide the set V^^^ into three subgroups: Vi^^^^ = V^^i, V^^^^^ = ^iv2 
and V^f^^3 = Vi^^ - V^^^^ - Vi^^^- Similarly, we divide the set V^^^. into three subgroups: 
"^UkI = Hvi^ HnK2 = ^yP3 and V^n,,3 = ^luj, - ^n,,! - ^UkI- Clearly, the subgroups 
VxQj^j, 1 < i < 3 are disjoint, so are true for subgroups VyQ^^j, I < j < 3. 

Denote j.^, ^.^ and ^.^ as the instantaneous sets of S^, V*.„ and when a transmitting 
node Xk is being explored with i = C{xk) in iteration + 1 of our algorithm. We divide the 
set into four subgroups as follows. Let S'^^^^^i = {{x,y) E Sl^^^,N'{x) E VtK,M{y) E 

^k}^ ^;.„2 = {{x.y) e ^;.,,Ar(x) E nK.U{y) e n^}, ^;.„3 = {{x,y) e si,^.U{x) e 

n^K,M{y) E fi^} and £l,^^, = {{x,y) E 4.,,Ar(x) E n\M{y) E ^k}. Divide V^,., and 
^iu,x^ respectively into four subgroups Vl^^^^^ and V*„^^^. ^ corresponding to ^^^^ 1 < j < 4. 
Clearly \£l^J + \£l^J + \£l^J + \£l,^J = K and rank(r(£;,J) = K for any ever 
explored by MDFS in iteration A' + 1 of our algorithm based on Lemma [3l 

In the following of this section, we prove a sequence of lemmas before we finally prove 
Theorem [2l Unless otherwise stated, we assume that we are in the last iteration X + 1 of our 
algorithm. 

Lemma 7: Let Xk be a transmitting node that has been explored by MDFS in iteration K +1 
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of our algorithm with C{xk) = i- When Xk is explored, we have Vflp^^ C ^.^ ^ + V*„ ^.^ 2 ^^'^ 
rank(r(VL,,^.,, + V:,,,^,^ + V;,,,J) (2) 

For any y G we have 

rank(r(Vi,,^,,i + VL,,^.,2 + x^, Vj„,,^ + 2/)) (4) 

= rank(r(V*„^^^^i + Vl.^^^^2-.^lu,xk + v)) = + l"^'L,Xfc,2l (5) 

Proof: Based on Lemma [6l all J\f{xj) with G V'xspan will finally be explored by MDFS 
in iteration K + 1 if Xk has ever been explored. Since for any x G V^^ 3 + V*„^^ 4, A/'(2;) 
is not explored by MDFS, so we have V^^^^^ C V^^ ,^^ 1 + V^^ ,,^ 2- By definition of V^|p„„, 
it's easy to conclude that Q-© hold. By definition Mdi) is not explored by MDFS for 
any y G V^j^^g. Based on Lemma [6l we have rank(r(V^„^^^ + Xfc, V^„^^^ + y)) = K. Given 
rank(T(Vl.„ ,,.^, V^„ ,,.^ + y)) = K and the fact that V^,*=p„„ C V^^^^ is the unique set satisfying 
rank(r(V;„,,^ + x,,Vl^^J) =rank(r(V;„,,^, V^„,,J) = K, we conclude rank(r(V;„,,^,i + 

Lemma 8: For any with C{xk) = i explored by MDFS in iteration + 1 of our algorithm. 

Proof: By definition all the super nodes where the nodes in V^-p^, V^j,^, Vy-pi and V*^3 
belong to are labeled unexplored when MDFS returns, so during the running time of MDFS in 
iteration K + 1, the nodes in V^-p^, ^xVA' ^Ivi ^^'^ ^Ivs never explored and V*.^3 C V*.„ j.^, 
^ivA ^ yL,x^^ ^iri ^ '^iu,xk '^yvz ^ ^iu,^u ^Iw^y^ hold for any Xk with £(2;^) = i explored 
by MDFS. Any A/'(x) (or M{y)) with transmitting node x (or receiving node y) used by but 
not by £p must be explored by MDFS, which means that for any Xk with = i explored by 

MDFS V*p34 = V*p3 + V*^4 is the complete set of transmitting nodes in ^.^ so that for each 
transmitting node x in this set, Af{x) is unexplored and Vy-p^^ = Vypi + V^^g is the complete set 
of receiving nodes in V^^ ^.^ so that for each receiving node y in this set, ^{y) is unexplored. So 

l^u,Xfe,ll + I'^u.Xfc.sl = I^T'll + I'^PsI = Ki + Ki-i and \£u,Xk,z\^\^u,XkA\ = I^T'sI + l'^P4l = ^i3 + ^i4- 

l^n,xfe,il ~ \^u,XkA\ ^ l^-pil ~ I^P4l = ~ -^«4 ^ Straightforward result by subtracting one 
from the other. ■ 
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Lemma 9: rank(T(VL,x„i + VL,x„2, "^yu.^uA + "^^",=^^,3)) > |VL,x„il - |VL,x„4l =Ki- 
for Xk being explored by MDFS. 

Proof: If i^j — iCj4 < 0, the statement is obviously true. Assume Ki — > 0. We know 

rank(r(£^*^2.^J) = K = + \Sl^^^^^2\ + \^u,xk,3\ + l'^M,Xfc,4l = rank(r(V*„^^^^i + V*„^2.^^2 + 

V* _I_V« v« -l-V* -l-V* -l-V* "l") TfrankfTfV* -l-V* V* -I- 

^yu,^„3)) < \^L,x„i\-\^L,x,aI we will have rank(T(V;„^^^^i+VL,x„2+VL,x„3+VL,x„4' 

Let Xfc be some transmitting node explored by MDFS in iteration K + 1 of our algorithm. 
Let y' G .j.^ 2 + ^yu,XkA- By definition, N{y') is explored by MDFS sometime. If y' G V*p, 
then y' must have been deleted from in a type 3 forwarding move of MDFS sometime in 
iteration K + 1 of our algorithm given that ^{y') is explored and have already been added 
back to if this type 3 forwarding move on y' happens before the current exploration of Xk 

otherwise it won't appear in 3,^. If y' ^ V*^, then y' must have been added to V*„^^ in a type 
1 forwarding move by MDFS before the current exploration to Xk otherwise y' won't appear in 
Vyu^xk when Xk is explored. It means for each y' G V*u ^.^ 2 + ^^n,xfe,4' v' either added to V*„^^ 
before the current exploration of Xi. or is deleted from V* , ^ after the current exploration of x^. 

Lemma 10: If rank(T(V;,,,^,i+VL,,^,2, V^«,x„i+V^„,,„3)) > |VL,x„i|-|VL,,„4l = ^.-^^4 
for Xk being explored by MDFS, then there exists some nonempty set Vy2A,Xk ^ '^lu,Xk,2+^lu,XkA' 
such that for any y" G V,24,.„ T(y^u,x,,i + Ku,x,,2,y") = E,.ev;' + ^xu,x,,2,y^) ior 

some = Vy2i,x, - y" + V^i with V^i C V^„^^^^i + V^^^^^^g. Let ^„24,x, ^ E^^^ .^ + ^ be 
the edge set with Vy2i,xk being their receiving nodes. 

Proof: We know that rank(T(V;„,,^,i + V;„,,^,2, Vl^,x„i + V^„,,„2 + V^„,,^,3 + Vi^^^^J) = 

\^xu,x,,,l\~^\^xu,Xk,2\- ''^^^^{'^{'^xu,x^,l~^^xu,x,,,2j^yu,x^,l~^^yu,Xk,3)) ^ I "^'xtt.Xfc,! I ~ I l^xw.Xj.^ I ' thcrC 

must exist some y' G V^„,,^_2 + V^„,,„4' such that T(V;,,,^,i + V;,,,^_2, ^0 = E,,ev; ^(^^,.,,1 + 
^L,x„2^yd for some C V^„^^^^i + V^„^^^,2 + V^„,^,,3 + V;„^^^^4 - y'. Let Vy24,x, = {y' + 
K} n {Ku,x,,2 + V^„,.„4}- Obviously Vy24,x, ^ and E,,ev,,4,., + VL,.,,2,2/.) = 

E,,ev;-v,.4,., nVL,x,,i + ^L,x,,2, yd- For any G V,24,.„ we have T(V:„,,^,i + V;,,,^,2, = 

VL,.„2, 2/0 with = V; - V,24,., C V;„,,,^,i + V;„,,,^„3 and V; = V,24,., - ?/" + V^. ■ 
Assume rank(r(V;„^^^_i + V;„,^^_2,V^„,x„i + V^„,^„3)) > \^L,x„i\-\^L,x,a\ =^i-^i4-Let 
y" £ l^?y24,a;fc bc the last ouc in Vy24:,xk being added to the set V*„^^ before the current exploration 
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of Xk or the first one in Vy24:,xk being deleted from the set ^.^ after the current exploration 
of Xk- By definition, 

+ VL,x,,2, y") = E T{VL,x,,i + Kn,.„2, y^) (6) 

for some = V,24,., - y" + C V;„^^^ - y" with C V;„,,^,i + V^„,,^,3. Let x" be the 
corresponding transmitting node such that (x",y") E £u24,x,,,- 

Lemma 11: Assume rank(T(V;„^^^_i + V;„,^^_2, + ^^«,x„3)) > WL,x^,i\ - WL^x^A = 

Ki — Ki4^ and let {x",y") E £u24,xk be defined above. Then just after adding y" or just before 
deleting y", we have 

^("l^xu.x",! + '^lu,x",2 + ^"i y") = ^C^'xu.x",! + '^xu,x",2 + (7) 

y^<^Vy 

for the same V^' as in Equation And when x" is explored just before adding y" along 
edge (x", y") or just after deleting y" along edge (x", y"), we have rank(T(V*„ 3,,,^ + V*„^„ 2 + 

Proof: Note that in equation (|7]), since y" E Vy2A,x^ with (x", j/") G is the last one in 
the set Vy24,Xk being added to the set V*„ ^.^ before the current exploration of x^ or the first one 
in the set Vy24,xk being deleted from the set V*„^^ after the current exploration of Xk, Vy24,xk~y" 
is not changed and is also subset of V*„^//. So in equation (|7]), Vy is the same as in © given 
that Vy24,xfe - y" is not changed and V^'i C ^ + V*„_^^^3 = V^p^ + V^^g is not changed, 

but the set V^ ,,^, 1 + V^^ ,,^ 2 are changing to V^.„^„ ^ + 2- When MDFS proceeds from 
the point of just after adding (x", y") to £u to the point of exploring Xk or from the point of 
exploring x^ to the point of just before deleting (x",y") from S^, only those three forwarding 
moves in definition of MDFS are allowed. It is sufficient for us to show that any forwarding 
moves of MDFS or backtracking of these moves doesn't change the relationship in equation Q 
so that equation ([7]) holds. 

Let's first consider forwarding moves of MDFS. A type 1 forwarding move along edge (x, y) 
would change V;„ ,, to V^^^^ + x. Since V^,p„„ C V^ ,, i + V^^ ,, 2 (based on Lemma |7]), the vector 
T(x, 3,) is a linear combination of row vectors in 7'(V*„^^ + V*„ 3, 2, 3,), so the relationship 
in ^ still holds while V^^^x,! + ^lu,x,2 changes to Vlu,x,i + ^xm,x,2 + a; in a type 1 forwarding 
move. In a type 2 forwarding move, V^ ,, i + V;„ ,, 2 changes to VL,x,i + ">^L,x,2 + a;-a;i for some 
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e V^span- Again we have Vf,p„„ C V^ ,, i + V^^ ,^ 2 (based on Lemma |7]). It is easy to see that 
the relationship in ^ still holds when V*.„^ ^ + V*„^ 2 changes to V*u 3, ^ + 3, 2 + x — Xj. In a 
type 3 forwarding move, some y G V*„^\Vj/24,Xfc is deleted from V*„^ which obviously doesn't 
affect the relationship in Q. Now let's consider backtracking moves of MDFS. Let equation Q 
hold after a type 3 forwarding move of MDFS along an edge (x, y) for some y ^ Vy2A,Xk- After 
the type 3 forwarding move along (x, y), MDFS will explore x when it explores J\f(x). Again we 
have V;^span ^ ^xu,x,i + ^L,x,2' equation © should hold before the type 3 forwarding move 
when V*„^ i + V*„^ 2 was Vlu,x,i + ^lu,x,2 + ^- A proceeding type 2 forwarding move before the 
current exploration of x means VL,x,i + ^L,x,2 was V^,:,. 1 + VL,x,2 + x - x' with x G Vflp„„ 
before the move. Let equation Q hold after a type 2 forwarding move of MDFS along a path 
P^'^^ with V;„^^^i + V;„^^^2- Again we have V^,p^„ C V^^^^^i + VL,x,2 (based on Lemma |7]), so 
equation Q holds with addition of row for x, which means that equation Q holds before the 
type 2 forwarding move when Vj^„ 1 + V^^^^^a was V^^^^i + Vl„ 2 + a; - x'. If the proceeding 
move is a type 1 forwarding move, it means that V^„^^ + V*„^ 2 was V*„^ ^ + 3, 2 - x, it is 
obvious that equation © holds before the type 1 forwarding move with V*„^ ^ + 3. 2 — x if 
it holds after the move with 3, ^ + V*„^ 2- From above discussion, we conclude that equation 
© holds. 

We know that after adding {x",y") or before deleting {x",y"), rank(T(V*„ ^ + V*^^/, 2 + 
^"1 ^lu,x"+y")) = \^xu,x",i\ + \^xu,x",2\ + '^- Given that equation ^ holds, we have rank(T(V*„ .j.,, 1+ 

'^xu,x",2 + '^lu,x")) = + I "^'x«,a;",2 1 + 1- * 

Lemma 12: For any x^ that has been explored by MDFS in iteration K +1 of our algorithm, 

rank(T(V*„ .j^^^ + VJ^u,Xk,2^^yu,Xk,l + "'^yu.Xfc.s)) ~ l"^'xn,Xfe,ll ~ l"^'a:u,Xfe,4l = ~ > 0. 

Proof: We first prove Ki — Kn > 0. Assume Ki — < 0. Then we must have 

rank(r(V:„,,^.,i + V;„,,^„2, V^„,,^,i + V^^^^J) > {V^^^^J - \V,^^,J = lU - K,,. Based on 
Lemma M and M if rank(r(V;„,,^,i + V;„,,^,2, ^lu,x,,i + ^U^^.^z)) > \^L,xJ ' l^kxj = 
Ki - Kii, then we will have rank(r(Vl^^„^i + Vi^y,^^ + x", V*„^^„)) = WL,x'\i\ + |VL,a;",2l + 1 
when x" is explored just before adding y" along edge (x", ?/") or just after deleting y" along 
edge {x",y") for x", defined as in Lemma [TT] But based on Lemma |71 we should have 
rank(r(V^„ ,^„ i + V^„^^„ 2 + x", V*„^^„)) = \^L,x",i\ + \^xu,x",2l which constitutes a contradiction. 
So we must have — Ki^ > 0. 

Assume rank(T(V;„^^^_i + V;„^^^,2, ^lu,x„i + Ku,x„3)) > \^xu,x„i\ - \^L,x„i\ = - > 
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0. Using a similar argument as above, we would arrive at a contradiction which means the 
assumption doesn't hold. So we must have rank(T(Vl.„^^^_i + Vi^^^^^2,Hu,x,,i + ^u,x^,3)) < 
l^xM,a:fe,il~I^L,a:fc,4l = ^i-Ki^. Now togcthcr with LemmajH wc concludc that rank(T( V*.„ ^.^ 1 + 

'^xu,Xk,2^^yu,Xk,l ~^ '^yu,Xk,3)) ~ \^xu,Xk,l\ ~ \^xu,XkA\ ~ ~ -^«4- ' 

Lemma 13: For any Xk with C{xk) = i that has been explored by MDFS in iteration K + 1 
of our algorithm and any G V^f^^g, rank(T(VL,a.„i + VL,x„2 + Xk,Viu,x„i + Hu,x„3 + 

%■)) =r^^HT{VL,x„i + ^L,x„2 + ^k, v^,,,^,i + v;„,,^,3)) = \vi^,,J - \vi^,,j = K,-K, 



i4- 



Proof: First we prove rank(T(V;„^^^_i + V^^^^^^g + Xk,Viu,x„i + Hu,x„3)) = K., - Ki^. 
Based on Lemma H V^'^p„„ C V^^^^^^i + V*„^^^^2' so rank(T(V;„^^^_i + V^^^^^^s + ^^fc, V^^^^^^i + 

'^yu,Xk,3)) ~^^^^i'^O^xu,Xk,l + '^xu,Xk,2^^yu,Xk,l + "^'j/w.Xfc.s)) ~ I ~ l"^'xu,Xfc,4l = Ki — Ki4 

based on Lemma [T2l 

Second we prove rank(r(VL,x„i + ^L,x„2 + Xk^iu^x,,! + V^„,x„3 + %)) = |VL,x„il - 
\^lu,xj = ^ - ^4. Given rank(T(V^„^,^^i + V^^^,^^ + x^, V^„,,^,i + Vl^^,J) = \Vl^^,J - 
|VL,x„4l = Ki- Kii, rank(T(VL^^^^i + VL,x„2 + Xk,Viu,x„i + Hu,x„3 + Vj)) equals either 
~ l"^'xu,Xfc,4l ~ l"^'xu,Xfc,4l + 1- Assume 

rank(T(V;,,,^,i + VL,,^,^ + ^k, V;,,,^,i + V;,,,^,3 + y,)) = \VL^,J - \V,,^,J + 1 (8) 
From Lemma |71 we have 

rank(T(V;„,,^,i + VL,.„2 + ^fc, K,,x, + %)) (9) 
= rmk{T{V;,,^,^^, + VL^,^^2,Ku,x,+y,)) = + |V;.,x„2l (10) 

and 

T{xk,Vl,^,^ + y,) (11) 
^ T(x',V;„,,^+%) (12) 

From dH]) and (fT2l) . we have 

rank(r(V*u 1 + .j,^ 2; "^'^n,^^,! + "^'yu.xfc.s + %■)) = ~ l"^'L,xfc,4l + 1 (13) 

From (flOl) and (fT3] ). we have 

r(Vk,.„i + ^kx„2,%) (14) 

I] r(vL,x„i + vL,x„2,?/i3)+ 5^ T{vL^,^^, + vL,x„2,y24) (i5) 

J/13eVyl3,a;j. J/24eVy24,3:j. 
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for some Vyu,^, C V^^^^^^i + V^„,^^ 3 and some nonempty Vy24,x, C V;„^^^^2 + V;«,x,,4- From 
(fT2l) and (fTSl), we have 



r(VL,x„i + VL,,,,2 + a;fc,%) (16) 

for the same Vyrs,xk Vy24,xk in (fT5l) . 

Let G Vy2i,xk l^st orie in Vy2i,xk being added to the set V^^^,^ before the current 

exploration of Xk or the first one in Vy24,x^, being deleted from the set after the current 

exploration of Xk and {x",y") E ^u^x^- Then using a similar argument as that in Lemma [TTl 
we have rank(T(V;„^^„^i + V*„^^„^2 + yUx" + Vj)) = + |VL,x",2l + 1 when x" is 

explored just before adding y" along edge {x",y") or just after deleting y" along edge {x",y"), 
but it's a contradiction with Lemma |71 So it must be rank(T(V*„^^ ^ + V*^ ^,^ 2 + ^fc, '^lu,xk,i + 

Lemma 74; rank(r^^) =rank(r(V:f,^, V^^^)) = - K,^,l < t < L - 1. 

Proof: If = 0, i.e., the cut and layer cut i has no intersection, then rank(Tj^^.) = 
rank(T(V;f^^, V^j^^ )) = Ki - Ki^ = holds. Next assume that S'^^ ^ 0. Lemma [T3] savs for 
any Xk with C{xk) = i that has been explored by MDFS in iteration A' + 1 of our algorithm 

and any yj G V^j^^g, rank(r(Vl^^^_^i + V^„^^^.^2 + Xk,Vi^,x„i + Ku,x„3 + Vj)) = Ki - K^^ and 
rank(r(VL,x^.^i + V;„^^^^2 + a^fc, V^«,x„i + ">^»,3)) = Ki-K^. LemmaHsays for any e V^j^^g, 
T{xk, Vy,^,^ + yj) = E.'gvS,.„cvi„,,^,,+Vi„,,^,, T{x', V^„,,^ + Based on these two Lemmas, 
it's easy to conclude that 

rank(T(VL,,^,i + VL,.,,2 + ^k, V;,,,^,i + V;„,,^,3 + Vjn^s)) (18) 
= rank(T(VL,,^,i + VL,,,.,2 + x,, V^^J) = - K,, (19) 
= rank(T(V:„,,^,i + VL,,^,2, V;^^,)) = K, - K,, (20) 



and 



nxk,Vl^,)= E ^(^'''^'^n,) (21) 



xu.Xu ,2 



Equations (1191 ) (1201) and (1211) hold for any x^ that has been explored by MDFS in iteration K + 1 
of our algorithm. 
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Let < q < Q he the gth transmitting nodes in layer i that has been explored by MDFS 
in iteration K + 1 in our algorithm. Note that since some transmitting nodes may be explored 
more than once, may not be distinct but Q is finite. We claim that 

rank(T(V;^i + T^p2 + E KnJ) = rank(T(K^i + V^^.,, V^^J) = - K,, (22) 

k=l 

for 1 < g < g. When is explored, V;„^i ^ = V^^^ and V*„^i 2 = ^xP2- From ^ we 
have 

rank(T(V:pi + V,^, + x\ V;^,)) = rank(T(V:pi + V^^^, V;^,J) = - lU, (23) 
When is explored, 

'^xu,x^,l + "^'a;n,x2,2 — '^xu,x^,l + "^'a:u,a;l,2 + ^ (24) 

and from (US]) ^ 

rank(r(VL,,.,i + VL,,.,^ + V;^,,)) = rank(T(VL,,.,, + VL,,.,^, V;^,)) = i^. - (25) 

From ^ to (|25]), we conclude that rank(r(V^pi + V*p2 + + ^^ V^j^k-)) = " ^i4- Use 
induction on a;"^, 1 < g < Q, just as we did for x^, we conclude that (l22l) holds for any 
g,l < g < Q. We know that V*,pi + V^pa + Efc=i a^'' = "l^xn^' ^o ^hen g = Q equation ([221) 
means that rank(r^^) =rank(T(V;f,^, V^^^^J) = - K,,,l < z < L - 1. ■ 
Lemma 15: 

L-1 L-1 
1=1 i=l 

Proof: By definition, Yl!i.=i the total times the paths in V cross the cut Qk. Since 

each path in V crosses Hk at least once and it's possible that some of these paths may cross 
Qk more than once, we have J2i=i — ^- For VP G V, let kp be the times P goes from 
VLk to Vfj^ and k'p be the times P goes from to VIk- Given 5* G VIk and D G ^2^, it 
must be true that kp — k'p = 1 and J2pev(^P " ^'p) ~ 1^1 ~ definition, the times that 

VP G P goes from to is counted in Ei=L^ and the times that VP G P goes from 
Vt%^ to fii^ is counted in Ya=i ^a, i.e., J^pevi^P ~ ^p) = J2i=i - Ya=i ^i^- So we have 

Et"i' - EtY = i^. ■ 

Based on Lemma \T4[ we have 

L-1 L-1 L-1 L-1 

mnk{T{Sn,)) = ^ rank(T^^) = ^(K, - K,,) = ^ i^, - ^ (27) 

i=l i=l i=l i=l 
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In Lemma [T5l we show that J2i=i ~ J2i=i — ^ which means rank(T(£'f7^,)) = K and 
this concludes our proof for Theorem [21 

■ 

Theorem [T] proves that our algorithm terminates in finite time. Theorem |2] proves that our 
algorithm finds C linearly independent S — D paths (V) where C is the unicast capacity of the 
underlying deterministic relay network. Lemma |4] shows that these C paths in V correspond to 
a capacity-achieving transmission scheme. They consist the complete proof of correctness for 
our algorithm for finding the unicast capacity of any linear layered deterministic wireless relay 
network. 

An arbitrary deterministic relay network Q can be unfolded over time to create a layered 
deterministic network Ql through time-expansion technique [1][2]. The transmission scheme in 
Ql identified by our algorithm corresponds to some equivalent transmission scheme in Q maybe 
time-variant. It means that our algorithm works for finding the unicast capacity of an arbitrary 
linear deterministic relay network. 

V. Conclusions 

The deterministic channel model for wireless relay networks has been a useful model for 
studying the capacity and capacity-approaching transmission schemes for underlying networks. 
In this paper, we proposed a fast algorithm for finding the unicast capacity of any given linear 
deterministic wireless relay network. Our algorithm finds the maximum number of linearly 
independent paths for the deterministic relay network and these paths correspond to a capacity- 
achieving transmission scheme. The essential component of our algorithm is a modified depth- 
first search algorithm developed for linear deterministic wireless relay networks. The proof of 
correctness for the algorithm is given which guarantees that our algorithm works in universal 
cases. Compared with previous results on solving the same problem, our algorithm prevails with 
significantly lower computational complexity. Moreover, the development of the modified depth- 
first search algorithm is based on a very intuitive idea, that is, to build up the path by adding 
edges while avoiding the linear dependency by using rank check at the same time. 
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